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SUMMARY 


This  report  presents  some  refinements  to  the  wing  weight 
estimation  methods  outlined  in  RAND  Report  R-100.(1)  Since 
Report  R-100  was  published,  the  weight  equations  have  been 
applied  to  a  number  of  modern  airplanes  to  determine  certain 
empirical  constants  and  to  observe  the  relative  size  of  the 
various  terms  in  the  equations.  Some  of  these  terms  were 
found  to  be  relatively  small.  Therefore,  it  was  possible 
to  shorten  the  final  equations  with  little  loss  of  accuracy. 

This  recent  study  also  showed  that  the  original  assumptions 
regarding  the  distribution  of  dead  weight  in  the  wing  needed 
revision.  A  more  accurate  method  of  accounting  for  the 
effects  of  dead-weight  loads  is  outlined  in  this  report. 


For  references,  see  p.39. 


SYMBOLS 


A  ~  area 

b  ~  wing  span  (aerodynamic),  in. 

bs  —  wing  span  (structural)  (=  6/cos  A),  in. 

Cc  —  intercept  value  for  t/l o  for  type  of  construction 
used  in  flanges  of  the  structural  box 

Cr  =  intercept  value  for  t/h  for  type  of  construction  used 
in  ribs 

Cs  =  intercept  value  for  t/h  for  type  of  construction  used 
in  the  shear  material  of  the  structural  box 

C  —  mean  aerodynamic  chord,  in. 

C;  —  integrated  effective  moment  coefficient 

Cs  =  structural  mean  aerodynamic  chord  (=C/ cos  A),  in. 

d=  ratio  of  thickness  taper  ratio  to  planform  taper 
ratio  ( =  m/A) 

Fa  =  equivalent  axial  allowable  stress  for  bendi  n  g ,  1  b/sq  in  . 

Far  =  allowable  compressive  stress  in  rib  flanges,  lb/sqin. 

h  —  maximum  wing  depth  at  any  span-wise  station,  in. 

J n  =  weight  of  leading  and  trailing  edge  structure  (in¬ 
cluding  flaps  and  ailerons)  divided  by  the  total  wing 
area 

Jn  =  inertia-load  correction  factor  (as  used  in  Report 
R-100) 

Jnb  =  inertia-load  correction  factor  for  bending  loads 

Jns  =  inertia-load  correction  factor  for  shear  loads 

Js  —  shear  reduction  factor  (as  used  in  Report  R-100) 

Jt  =  integrated  shear  reduction  factor 

Kg  =  width  of  wing  structural  box  expressed  as  a  fraction 
of  the  chord 

K i  =  rib  spacing  expressed  as  a  fraction  of  chord  (=L/C) 

Ky,  =  maximum  wing  thickness,  expressed  as  a  fraction  of 
chord 

kL  =  ratio  of  effective  column  length  to  actual  length 

(=L0/L) 
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k  ■ 

Rlbd 

*i . 

kita 

ki,d 

krf 


span-distribution  factor  (as  used  in  Report  R-100) 

sp an - d i s t ri bu t ion  factor  for  airloads  (also  equal 
to  fe  *  ) 

span- dis tribu tion  factor  for  inertia  loads  (also 
equal  to  kitd  ) 

wing  effective  depth  factor 

bending  integration  factor  (as  used  in  Report  R-100) 
bending  integration  factor  for  airloads 
bending  integration  factor  for  inertia  loads 
shear  integration  factor  (as  used  in  Report  R-100) 
shear  integration  factor  for  airloads 
shear  integration  factor  for  inertia  loads 

ratio  of  rear  spar  depth  to  airfoi)  maximum  depth 

(=  hrf/h) 

ratio  of  average  spar  depth  to  airfoil  maximum  depth 
(  =  ^aveA) 

nonoptimum  factor  for  wing  primary  structure 

nonoptimum  factor  which  includes  primary  and  secondary 
structure 


M  =  bending  moment,  in. -lb 
Ma  =  bending  moment  due  to  airloads,  in. -lb 
Mj  —  bending  moment  due  to  inertia  loads,  in. -lb 
m  =  thickness  taper  ratio  (=  hT/hR) 
mar  =  ratio  of  allowable  rib  flange  stress  to  equivalent 
axial  stress  (=  WM 

ms  =  ratio  of  effective  shear  stress  to  equivalent  axial 
stress  Fs0/Fa) 
r\j  —  ultimate  load  factor 
P  —  load,  lb 

qv  —  dynamic  pressure,  lb/sq  ft 
5  =  wing  area,  sq  ft 

tr  —  minimum  average  thickness  for  ribs,  in. 

Va  =  shear  due  to  airloads,  lb 
Vj  =  shear  due  to  inertia  loads,  lb 

=  value  of  shear  measured  at  the  root  section,  lb 
(fg  =  airplane  gross  weight,  lb 
If.  =  weight  of  wing,  lb 
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=  weight  of  wing  and  wing  contents,  lb 

)]  =  nondimensional  span-wise  distance  measured  from  the 
root  section 

A=  planform  taper  ratio  (=  Cj/Cr) 
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I.  INTRODUCTION 


RAND  Report  R-100  outlined  procedures  for  applying  optimum 
design  methods  to  structural  weight  analysis  of  the  wing. 
Subsequent  to  the  publication  of  Report  R-100,  the  wing 
weight  equations  were  applied  to  a  number  of  modern  air¬ 
planes  to  check  the  values  of  calculated  weight  against  the 
actual  weight.  This  preliminary  study  showed  the  need  for 
revision  of  some  of  the  original  simplifying  assumptions 
regarding  the  distribution  of  dead  weight  within  the  wing. 

The  procedure  for  estimating  wing  weights  may  generally 
be  analyzed  in  two  parts: 

1.  The  estimation  of  the  loads  applied  to  the  structure, 
based  on  a  knowledge  of  the  geometric  properties  of 
the  wing  and  of  certain  loading  parameters 

2.  The  estimation  of  the  weight  of  structure  necessary 
to  carry  these  loads. 

The  accuracy  with  which  weights  can  be  estimated  depends 
diiecciy  on  the  accuracy  with  which  the  designing  loads  are 
estimated.  The  revision  of  the  integration  factors  represents 
a  means  of  improving  the  results  of  the  weight  equations  by 
use  of  a  more  accurate  method  of  determining  the  wing  design 
loads . 

For  Report  R-100,  the  airloads  on  the  wing  structure  were 
determined  by  use  of  Weissinger’s  method. <2)  These  loads  in 
turn  were  used  to  calculate  integration  factors  as  functions 
of  planform  taper  ratio,  sweepback,  and  thickness  taper 
ratio.  The  integration  factors  represent,  in  nondimensional 
form,  the  volume  of  material  necessary  to  carry  the  airload 
shear  or  bending  moment  for  a  given  wing  configuration.  How¬ 
ever,  to  obtain  the  net  shear  and  bending  moments,  the 
inertia  loads  must  be  taken  into  account.  This  may  be  done 
in  either  of  two  ways:  (1)  by  subtracting  inertia  loads  from 
the  airloads  and  working  with  net  loads  only,  or  (?)  by 
calculating  integration  factors  separately  for  airloads  and 
for  inertia  loads  and  using  these  factors  to  calculate  the 
net  result.  The  second  of  these  two  methods  is  used  here 
because  of  its  greater  convenience. 


For  references,  see  p.39. 
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Another  revision  made  in  the  original  equations  is  a 
change  in  the  terms  representing  rib  weight.  These  changes 
have  been  included  in  order  to  make  the  equations  more  accu¬ 
rate.  Derivation  of  the  new  terms  is  given  in  Appendix  C. 

When  the  original  equations  of  Deport  R- 100  were  tried 
out  on  wings  corresponding  to  the  configurations  and  loadings 
of  several  modern  wings,  it  was  found  that  some  of  the  terms 
could  be  neglected  with  no  appreciable  error.  Therefore  it 
was  possible  to  write  the  wing  weight  equation  in  a  shortened, 
more  convenient  form.  This  shortened  form  of  the  equation  is 
given  in  Section  III. 
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II.  CHANGES  IN  THE  ORIGINAL  WING  WEIGHT 
EQUATIONS 


The  expression  for  the  volume  of  the  optimum  wing  struc¬ 
tural  box  material  is  given  in  Report  R-100  as 

(  ,  ,  -  KBCiqv  - 

Vol  =  2 C,k*K£  +  CckLKLKB  +  Cr  ~~~  )  1 44CS  +  2  CS 

2  Fa  (J2  ms  3  16  ke  hB]  ms  J 
(See  Report  R-100,  Eq.89.) 

This  equation  comprises  the  following  expressions  for  the 
volume  of  the  different  structural  components: 

Shear  Material 

Vol  =  2C$  ksK^(144CS)  size  term  (la) 

+  —  '*r'nf  8  loading  term  (lb) 

2  Fa  \2  ••  *sy 


Bending  Material 
Vol  =  CckLKLKB{  144CS) 

+  1  JnnfWgb  ( 1  4  V2\ 

2  F.  \3  h,j 


size  t  e  rm 
loading  term 


(  lc) 
(Id) 


rib  flange  material 


One  part  of  Eq.(l)  found  to  need  revision  is  the  ex¬ 
pression  for  the  volume  of  rib  material.  This  change  is  made 
so  that  the  equation  describes  more  accurately  the  present 
practices  in  rib  construction. 

In  order  to  provide  a  minimum  gauge  restriction  on  the 
rib  shear  material  rather  than  the  minimum  proportion  re¬ 
striction  provided  by  Eq.(l),  the  volume  of  the  size  term 
for  rib  shear  material  is  rewritten  as 


V  1  =  *r0  KBKhS(lU) 

Kl 


(2) 


where  tr<>  =  minimum  average  thickness  of  rib  shear  material. 
The  derivations  for  the  new  rib-material  terms  are  presented 
in  Appendix  C. 

The  expression  for  the  volume  of  rib  flange  material, 
given  by  Eq.(lg),  is  ->termined  from  a  consideration  of  the 
amount  of  flange  material  necessary  to  resist  chordwise 
bending  moments  from  the  trailing  edge  structure.  Since  some 
of  the  parameters  in  Eq.(lg)  may  be  difficult  to  evaluate 
in  preliminary  design  stages,  the  volume  of  rib  flange  mate¬ 
rial  will  be  based  on  a  different  set  of  parameters. 

The  revised  expression  for  rib  flange  material  presented 
in  E  q . ( 3 )  gives  the  amount  of  material  necessary  to  carry 
the  axial  loads  associated  with  chordwise  bending  moments. 
In  this  case,  the  bending  moments  result  from  the  transfer 
of  normal  airloads  from  the  point  of  application  in  a  chord- 
wise  direction  to  the  shear  resisting  material  in  the  spars. 
The  expression  for  the  volume  of  flange  material  is 

_  Wg  (3) 

F.r 


where  Far  =  allowable  axial  stress  in  the  rib  flanges.  The 
derivation  of  Eq.(3)  is  given  in  Appendix  C. 

When  Eqs.(2)  and  (3)  are  substituted  for  the  corresponding 
rib  terms  ( le )  and  (Ig)  in  Eq.(I),  the  volume  of  wing  struc¬ 
tural  material  becomes 


Vol  =  [2CtklKjl  +  CckLKLKB  + 


\ 


trpWh) 

KlC  ) 

1  •V'/'t /T‘,  it 

2  f.  U2'--3* 


kb  (6/2) 

h h 


b  +  KgC 


Kb 


2  Khma 


1 

-+  — 


} 


(4) 
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where 


Bor 


The  term  (1/),  representing  the  loading  term  for  rib  shear 
material,  is  not  changed. 

The  optimum  weight  of  the  wing  structural  box  is  found 
by  multiplying  Eq.(4)  by  the  density  w.  Using  a  single  over¬ 
all  nonoptimum  factor  kx  and  adding  the  weight  of  leading 
and  trailing  edges  JltS,  the  total  wing  weight  is  written 


»'  = 


7  7  t  K  ftK  h 

.  2  £/  2  j_  /-»  l  1/  1/  -L  r0  n  n 


2  C,km  +  CckLKLK„  + 


KLC 


144CS 


X  1  Jn"fWE 
2  Fa 


(kis  Js  kih  kb  (6/2)\ 


2  ms 

+  KBC 


3  k . 


/  Kb 

+  4 

\2Kflm.ar 

ms) 

+  JLTS  (5) 


where  w  =  density  of  the  material 

k  _  actual  weight  of  structural  box 
x  optimum  weight  of  structural  box 

JltS  =  weight  of  leading  and  trailing  edge  structure 
including  flaps  and  ailerons. 
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III.  SHORTENED  FORM  OF  WING  WEIGHT  EQUATION 


The  terms  of  Eq.(5)  were  evaluated  for  several  wings  of 
contemporary  commercial  and  military  airplanes.  It  was  found 
that  the  terms  representing  rib  flange  material  and  the 
loading  term  for  the  rib  shear  material,  Eqs.(3)  and  (1/), 
could  be  omitted  with  no  appreciable  effect  on  the  final 
numerical  result. 

The  nonoptimum  factor  kx  is  defined  as  the  ratio  of  the 
actual  structural  box  weight  to  the  weight  of  the  optimum 
structural  box.  It  was  originally  intended  to  determine 
values  of  kx  by  comparing  calculated  values  of  Eq . ( 5 )  with 
actual  values  of  weight  for  the  airplanes  studied.  However, 
much  of  the  weight  data  available  at  this  time  was  not  item¬ 
ized  in  sufficient  detail  to  evaluate  Jir  with  the  required 
accuracy.  Therefore,  a  new  nonoptimum  factor  was  defined  and 
used  in  connection  with  the  statistical  study. 

To  avoid  introducing  any  additional  error  by  using  un¬ 
certain  estimates  for  values  of  Jlt,  the  new  nonoptimum 
factor  is  defined  so  as  to  include  the  weight  of  secondary 
structure  (leading  and  trailing  edge  structure). 

If  the  assumption  is  made  that  the  weight  of  the  secondary 
structure  is  proportional  to  the  weight  of  the  primary 
structure,  the  wing  weight  is  given  by 

^to  =  (kx  +  )  X  (  optimum  weight  of  structural  box).  (6) 
Let  the  new  nonoptimum  factor  be  defined  as 

k'x  =  kx  +  Ki  .  (7) 

Then  the  over-all  weight  becomes 

=  fe^X(optimum  weight  of  structural  box).  (8) 

In  this  form,  kx  could  be  readily  evaluated  and  was  found 
to  be  a  function  of  planform  taper  ratio,  as  predicted  in 
Report  R-100.  Due  to  the  classified  nature  of  the  weight 
data  involved,  actual  values  are  not  included  in  this  un¬ 
classified  report. 


I 

» 


Omitting  the  rib  weight  terms  discussed  previously,  the 
shortened  form  of  the  wing  weight  equation  is  now  written 


(■ 


2  Ctk*Kl  +  CckLKLKB  +  ImJ^JlUCS 

+  +  l  K  <6/2>V 

2  Fa  \2  l**s  3  *6fee  h«  ) 


(9) 


,  ,  ,  _  _ actual  wing  weight 

where  A*  —  optimum  weight  of  structural  box  * 
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IV.  INERTIA  LOAD  RELIEF 


As  shown  in  Report  R-100,  part  of  the  wing  structural 
material  is  proportional  to  the  size  of  the  wing  and  part 
is  proportional  to  the  loading.  The  integration  factors 
represent,  in  nondimensional  form,  the  volume  of  material 
which  is  proportional  to  loading.  These  factors  are  calcu¬ 
lated  separately  for  shear  material  and  bending  material. 

The  value  of  an  integration  factor  is  determined  only  by 
the  span-wise  distribution  of  shear  load  or  axial  flange 
load.  It  was  assumed  in  Report  R-100  that  the  dead-weight 
loads  had  the  same  distribution  as  the  airloads.  Therefore, 
the  distribution  of  net  loads,  when  plotted  nondimensionally, 
was  the  same  as  for  airloads  and  consequently  had  the  same 
value  of  integration  factor.  All  that  was  necessary  in  order 
to  include  the  effects  of  inertia  loads  was  to  correct  the 
value  of  loads  at  the  root  section.  This  was  done  in  Report 
R-100  by  correcting  the  shear  at  the  root  by  multiplying  by 
the  factor  J„,  where 

Jn  =  wJ-~  ^  (10) 

*1 

and  where  Hg  =  airplane  gross  weight 

y/wc  =  weight  of  wing  and  contents. 

(See  Eq.lOa,  Report  R-100.) 

Since  it  is  no  longer  assumed  that  the  inertia  loads  are 
distributed  in  the  same  manner  as  the  airloads,  the  simple 
correction  factor  shown  above  cannot  be  used.  A  study  of 
modern  airplanes  has  shown  that  the  difference  in  values  of 
airload  integration  factors  and  net  load  integration  factors 
may  be  of  considerable  magnitude  in  some  cases.  Therefore, 
a  method  was  derived  for  treating  the  airload  and  inertia¬ 
load  integration  factors  separately  and  combining  them  later 
to  obtain  net  values.  The  method  is  presented  below. 

The  integration  factors  are  obtained  from  integration  of 
the  loads  expressed  in  nondimensional  form.  For  shear,  the 
factor  is  obtained  from  an  integration  of  the  shear  loads. 
For  bending  material,  it  is  obtained  from  an  integration  of 
the  axial  flange  load  expressed  nondimensionally.  Since 
the  portion  of  the  material  considered  here  is  operating 
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at  constant  allowable  stresses,  the  integration  factor  also 
represents  the  amount  of  material  necessary  for  carrying 
the  shear  or  bending. 

The  airload  integration  factors  have  been  described  previ¬ 
ously  in  Report  R-100.  The  total  inertia-load  integration 
factors  may  be  composed  of  several  factors  which  are  obtained 
separately  for  concentrated  loads  and  for  distributed  loads 
and  which  are  then  combined.  A  procedure  for  the  calculation 
of  inertia-load  integration  factors  is  outlined  in  Section  V. 

After  the  airload  and  total  inertia-load  integration 
factors  have  been  obtained,  they  are  used  to  calculate  values 
of  the  inertia  relieving  load  factor  Jn .  This  factor  is  de¬ 
termined  separately  for  shear  and  bending  material.  As 
derived  in  Appendix  A,  the  expressions  for  Jn  are  found  to 
be  as  follows: 


Shear  Material , 


1 


VC 


(11) 


where  kisd  =  shear  integration  factor  for  inertia  loads  (total) 
kita  =  shear  integration  factor  for  airloads; 

Bending  Material , 

Jnb  =  1 


where  kbd  =  span-wise  load  distribution  factor  (  =  for 

inertia  loads 

kba  —  span-wise  load  distribution  factor  (  =  kisa )  for 
airloads 

kibd  ~  bending  integration  factor  for  inertia  loads 
kiba  =  bending  integration  factor  for  airloads. 

The  span  distribution  factor  kb  is  identical  with  the  shear 
integration  factor  k^.  Therefore,  values  calculated  for  kis 
may  be  used  in  Eq.(12). 

The  two  terms  ‘til  and  Jnb  ate  substituted  in  the  wing 
weight  equation  in  place  of  the  single  value  of  Jn  previously 
used. 

The  shear  reduction  factor  Jg  also  required  revision  in 
view  of  the  more  general  assumptions  made  regarding  inertia 
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load  distribution.  The  derivation  of  Appendix  B  gives  the 
expression  for  the  new  shear  reduction  factor  as 


j<  =  l  —  ^jbaSl.  ^nb 

s  3t  J 

J  Kg  *  _ 


.  ,  .  ,  tip  thickness 

where  *  =  thickness  taper  ratio  =  root  thickness 

ke  =  effective  depth  factor. 

The  corrected  wing  weight  equation  is  obtained  by  substi¬ 
tuting  the  new  values  of  Jn  and  Js  in  Eq.(5)  of  Section  II: 
r 

K  =  ^2C,A,altf  +  CckLKLKB  +  trft^^l44CS 

+  2  nf **  ["/i  Jn'kig.J',  .  2  ,  .  feba  (fc/2)\  , 

2  Fa  |\2  iJnb  l>‘  kg  hB)b 

+  +  — H  +  JtT5  .  (14) 

mar  J 

Substituting  the  new  values  of  J„  and  Js  in  the  shortened 
form  of  the  wing  weight  equation,  given  by  Eq.(9),  yields 
the  expression 


K  =  k'xw  2C,*Xa  +  CckLKLKB  + 


troHhKB 


144CS 


In  the  case  of  a  sweptback  wing,  the  values  of  structural 
span  and  structural  mean  aerodynamic  chord  are  used.  The 
structural  span  6,  is  defined  as 


bs  = 


where  A  =  sweepback  angle  measured  at  the  quarter-chord 
position. 

The  structural  mean  aerodynamic  chord  C*  is  defined  as 

Ct  —  C  cos  A . 


(17) 
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Equation  (17)  gives  exact  values  only  for  nontapered 
wings.  However,  except  for  low-aspect-ratio  highly  tapered 
designs,  these  expressions  give  sufficient  accuracy  when 
used  with  tapered  wings. 

In  order  to  write  Eq.(15)  in  a  more  general  form,  including 
the  effects  of  sweepback,  let 


P'=  ^2  C.Khh 


t.Kk  Kg 

+  CckLKLKB  + 


klcs 


?) 


e  'U- ^7“  3  “■  "‘^7  ~) 


(18) 

(19) 


Substituting  in  Eq.(15),  the  wing  weight  becomes 


V.  METHODS  FOR  ESTIMATING  INERTIA-LOAD 
INTEGRATION  FACTORS 

For  estimating  the  integration  factors,  the  inertia  loads 
are  broken  down  into  two  general  types — distributed  and  con¬ 
centrated.  The  distributed  loads  come  from  the  wing  structure 
and  from  miscellaneous  distributed  items  such  as  control 
systems.  The  concentrated  loads  originate  with  wing-mounted 
engines,  with  landing  gear,  and  sometimes  with  fuel  tanks.  It 
is  usually  more  convenient  to  treat  the  fuel  tanks  as  several 
concentrated  loads  rather  than  as  a  distributed  load,  since 
a  restricting  assumption  is  made  regarding  the  manner  in 
which  the  distributed  loads  vary. 

Integration  factors  are  calculated  separately  for  the 
different  loads  and  are  then  combined  into  a  single  factor. 
These  factors  are  calculated  separately  for  shear  loads  and 
for  bending  loads. 

CONCENTRATED  LOADS 

The  expression  for  the  shear  integration  factor  is  derived 
below  for  the  case  of  several  concentrated  loads  located  at 
different  span-wise  stations.  This  loading  condition  is 
shown  in  F ig. 1 . 

The  shear  curve  for  this  loading  is  shown  in  Fig. 2. 


Fig.  I — Concentrated  loading  condition 


Fig.  2 — Shear  diagram  *e»r 
concentrated  toads 


The  area  under  the  shear  curve  of  Fig. 2  is  given  by 
Area  =  Pifli  +  Pj  aj  +  P„  a„  . 


(21) 
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This  may  be  expressed  in  nondimensional  form  by  dividing 
by  the  value  of  shear  at  the  root  and  by  the  semispan  (see 
Appendix  F,  Report  R-100): 

Are.<  ,22) 

(ft+P»+P»)(V 2)  ' 

Since  the  allowable  shear  stress  is  a  constant  for  this 
portion  of  the  material,  Eq.(22)  also  represents  the  volume 
of  shear  material  expressed  in  nondimensional  form. 

The  integration  factor  for  shear  is  defined  as  the  ratio 
of  the  area  under  the  shear  curve  to  that  under  the  shear 
curve  for  the  basic  case.  Since  the  area  equals  one-half* 
for  the  basic  case,  the  shear  integration  factor  for  concen¬ 
trated  loads  is  given  by 


(6/2) 


Pl  a !  +  P2  a  2  +  Pn  an 
Pi  +  P2+Pn 


The  integration  factor  for  bending  depends  on  the  distri¬ 
bution  of  the  axial  flange  load  Pa ,  where 

M 

Pa  =  7  (24) 

h 


and  where  M  =  bending  moment 

h  =  depth  of  the  wing. 

The  integration  factor  fejj 
follows: 


for  bending  is  defined  as 


Let  A0  =  area  under  the  curve  of  M/h  versus  distance  along 
the  semispan  for  the  given  loading  condition 

•4ft  =  area  under  the  curve  of  M/h  versus  distance  along 
the  semispan  for  the  basic  case. 


Then 


kib 


(see  Report  R-100,  page  47.) 

For  convenience,  these  areas  are  expressed  in  nondi¬ 
mensional  form.  Since  the  area  (nondimensional)  under  the 
curve  of  M/h  for  the  basic  case  equals  one-third,*  the 
bending  integration  factor  is  written 

ki  =  34q  (26) 


*  See  RAND  Report  R-100,  pp.46  end  47. 
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where  9  =  1  —  T),  where  7]  is  the  nondimens ional  distance  from 
the  root  to  the  span-wise  station  at  which  the 
concentrated  load  is  located 
_  tip  thickness 
1,1  —  root  thickness  * 

A  plot  of  Eq.(28)  is  shown  in  Fig. 3.  A  value  of  is  read 
from  the  curves  for  each  concentrated  load.  These  separate 
values  are  then  combined  into  a  single  integration  factor 
by  the  expression 


0  0.1  0.2  0.5  0.4  0.5  0.6  0.7  0.8  0.9  1.0 


Fig.  3 — Integration  factors  for  bending  material — concentrated  loads 
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where  ki6  =  value  read  from  the  curves  of  Fig. 3 
P  =  value  of  concentrated  load 

T\p  —  nond  intension  al  distance  from  the  location 
load  to  the  wing  root  section. 

Equation  (29)  is  a  special  case  of  Eq.(33),  which  is 
in  Appendix  D. 


of  the 
derived 


DISTRIBUTED  LOADS 


The  integration  factors  for  distributed  loads  are  derived 
on  the  basis  of  the  following  assumption — that  the  load  at 
any  span-wise  station  is  proportional  to  the  cross-sectional 
area  of  the  wing  at  that  station.*  Such  a  load  distribution 
would  result  from  the  exact  filling  of  the  enclosed  volume 
of  the  wing  by  a  material  of  constant  density. 


The  shear  at  any  station  is  proportional 
wing  volume  outboard  of  that  station.  An  alg 
for  the  shear  integration  factor  yields  the 

k  =  1  +  X  +  d\  +  3d\2 
l,d  2+X  +  dX+2dX2 


to  the  enclosed 
ebraic  solution 
express  ion 


(30) 


where  d  —  m/X.  A  plot  of  this  function  is  shown  in  Fig. 4. 


0  01  0  8  0  3  0.4  0.5  06  0.7  0.6  0.9  10 

Fig.  4 — Shear  integration  factor  for  weight  distributed  as  wing  volume 


'Reference  3  has  shown  that  the  span-wise  distribution  of  structure 
weight  for  tapered  wings  tended  to  vary  as  the  square  of  the  chord.  For 
a  wing  having  a  constant  airfoil  section  along  the  span,  such  a  distri¬ 
bution  would  also  vary  as  the  wing  cross-sectional  area. 
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Similarly,  the  solution  for  the  bending  integration  factor 
is  given  by 


k'bd  B 


aX/1 


T  U"a+fl  ^  a 


l+a\  m— 2n\—  2/1  a 


- +  a2 


6(  1— m)  \3  2 


— 


a  loge 


l  +  a\  1~X/1  a  a2  .  .  1  +  a' 

-  +  - - + - a3  +  a4  log.  - 

a  /  12  \4  3  2  a  j 


(31) 


where 


B  = 


mX  la—  2/nX+X  (1—  X)  (1  —  at) 

-  +  - -  -+*  - 

2  6  12 


m 

a  =  -  . 

1  ~  m 


x 

Fig.  5  —  Integration  factors  for  bending  material  for  inertia  loads 
distributed  as  wing  volume 


If  it  becomes  desirable  to  revise  the  assumptions  regard¬ 
ing  the  variation  of  distributed  inertia  loads  in  the  wing, 
integration  factors  can  be  derived  for  any  assumed  variation 
of  distributed  loads.  These  factors  would  then  be  used  in 
the  same  manner  as  those  derived  here. 
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TOTAL  COMBINED  VALUES 


After  values  of  the  integration  factors  have  been  calcu¬ 
lated  separately  for  concentrated  and  distributed  loads, 
these  values  must  be  combined  to  give  total  values  for 
inertia-load  integration  factors.  This  is  done  for  the  shear 
and  bending  factors  by  means  of  the  expressions  given  below. 

The  total  value  for  the  inertia-load  shear  integration 
factor  is  given  by 


fei  .  = 

l*a 


_  kiA^VR+ki 


V*>  +  ^ 


+  kt  VH 

1  in  n  r 


+  v» 


(32) 


where  k{t  =  shear  integration  factor  for  any  given  load. 

Vr  =  corresponding  value  of  shear  at  the  root. 

The  total  integration  factor  for  bending  material  is 
given  by 


* 


MB1+kiiaMR2  +  •••  +  kj  bnMK 

Mrx  +  +  *  *  •  + 


(33) 


where  =  bending  material  integration  factor  for  any  given 
load 

Mr  —  corresponding  value  of  bending  moment  at  the  root. 
The  derivation  of  these  expressions  is  given  in  Appendix  D. 
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VI.  CONCLUSIONS 


In  order  to  improve  the  accuracy  of  the  equation  for 
estimating  wing  weight  it  was  necessary  to  use  a  more  accu¬ 
rate  method  of  estimating  the  loads  on  the  wing  structure. 
This  required  a  more  precise  procedure  for  estimating  the 
effects  of  dead-weight  loads,  such  as  the  method  outlined 
in  this  report.  The  use  of  this  refined  method  is  more  im¬ 
portant  for  large  airplanes  having  wing  fuel  and  wing-mounted 
engines  than  for  airplanes  which  have  no  sources  of  large 
concentrated  loads  in  the  wings. 

The  shortened  form  of  the  weight  equation  (Eq.20)  using 
the  over-all  nonoptimum  factor  fe*  is  the  most  convenient  form 
to  use  in  the  estimation  of  wing  weights  on  the  basis  of 
currently  available  weight  data.  The  nonoptimum  factor  fc*  may 
be  determined  from  statistics  and  expressed  as  a  function  of 
planform  taper  ratio. 

Values  of  airload  and  inertia-load  integration  factors 
were  obtained  from  actual  bending  moment  curves  and  were 
checked  against  values  calculated  by  the  methods  outlined  in 
this  report.  The  values  were  in  good  agreement,  the  best 
correlation  being  obtained  for  the  larger  airplanes. 
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APPENDIX  A 


DERIVATION  OF  EXPRESSIONS  FOR  INERTIA-LOAD 
RELIEF  FACTORS 


In  the  derivation  of  integration  factors  for  the  shear 
and  bending  materials  (see  Report  R -100)  it  was  assumed  that 
the  dead  weight  of  wing  and  contents  had  the  same  distri¬ 
bution  as  the  airloads.  The  use  of  Jn  as  outlined  previously 
gives  the  correct  net  values  of  shear  and  bending  moment  at 
the  wing  root  section  but  does  not  account  for  the  actual 
distribution  of  dead  weight. 

An  investigation  of  the  distribution  of  dead-weight  shear 
and  bending  moments  for  actual  airplanes  has  shown  that  the 
use  of  Jn  alone  yields  somewhat  optimistic  values  for  the 
weights  of  shear  and  bending  material.  This  is  particularly 
true  for  large  airplanes  having  engines  mounted  in  the  wings. 

In  order  to  make  the  wing  weight  equations  more  accurate 
by  accounting  for  the  actual  distribution  of  dead  weight,  the 
following  derivation  is  presented. 

The  shear  and  bending  materials  discussed  here  are  those 
portions  only  which  vary  with  loading  and  which  operate  at 
constant  allowable  stresses.  The  portions  which  are  pro¬ 
portional  to  size  are  not  affected. 


SHEAR  MATERIAL 

Equation  (84)  of  Report  R-100  gives  the  expression  for 
the  volume  of  shear  material  as 

(Vol),  =  (34) 

where  =  integration  factor  derived  on  the  basis  of  air¬ 
loads 

At  —  area  of  shear  material  required  at  the  root 


section  and  given  by 

i  _  1.  nf  Js  Jn 
*  '  2  * 


(35) 


Substituting  Eq.(35)  in  Eq.(84),  the  volume  of  shear  ma¬ 
terial  becomes 


(Vol),  =  . 


4* 


(36) 


*o 
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The  use  of  Jn  gives  the  required  area  of  shear  material 
corrected  for  the  relieving  shear  from  dead  weight.  However, 
in  actual  airplane  wings,  the  value  of  shear  relief  given  by 
the  term  Jn  is  correct  only  at  the  root  section,  since  the 
dead-weight  shear  is  not  proportional  to  the  airload  shear 
at  all  span-wise  stations.  This  difference  can  be  accounted 
for  by  using  a  combination  of  two  integration  factors — one 
for  axrioads  and  one  for  the  relieving  dead-weight  loads. 

The  net  cross-sectional  area  of  shear  material  required 
at  any  station  is  given  by 

X.  =  5Z*  (37) 

Fso 

where  Va  =  shear  from  airloads 

Vj  —  shear  from  dead-weight  loads. 

The  volume  of  shear  material  is  found  by  integration  of 
Eq. ( 37  ) : 

(Vol),  =  -p  fVadT)  -  fVrfdrj  .  (38) 

^s0  «/°  rs0  */° 

The  subscript  R  designates  values  taken  at  the  root 
section,  and  the  value  of  Va  at  the  root  section  is  given 
by 

%  =  lhnfJs  •  (39) 

The  value  of  Vd  at  the  root  section  is  given  by 

=  2  W 'vc  nf  Js  ( 40 ) 

where  c  =  weight  of  wing  and  contents. 

Multiplying  and  dividing  the  first  and  second  terms  of 
Eq.(38)  by  Eqs.(39)  and  (40),  respectively,  the  volume  of 
shear  material  may  be  written 


(Vol), 


_  *g  nf  J$ 


Jo  drx  -  W»cnfJ*b  Cl  il 

Jo  Vv 


This  assumes  that  the  shear  from  airloads  is  greater  than 
the  dead-weight  shear  at  all  span-wise  stations  and  was  true 
for  the  most  critical  design  condition  for  all  of  the  air¬ 
planes  studied. 


Since  the  allowable  shear  stress  is  constant,  the  cross- 
sectional  area  of  shear  material  is  directly  proportional  to 
the  shear  load.  Therefore, 


(42) 


where  A,a  —  area  of  shear  material  which  would  be  required 
to  carry  airloads  only.  Also, 


Vd_  =  A* 

%  V 


(43) 


where  ASd  =  area  of  shear  material  which  would  be  required 
to  carry  dead-weight  loads  only.  Substituting  Eqs.(42)  and 
(43)  into  Eq . ( 41 )  and  factoring, 


(Vol)s 


nfJs 


2  Fs 


(44) 


Since  the  integration  factors  are  evaluated  here  for  both 
airloads  and  inertia  loads,  it  is  necessary  to  use  subscripts 
a  (airload)  and  d  (dead-weight  load)  with  the  integratibn 
factors.  The  value  of  ^iJ0  of  this  report  corresponds  to 
kis  of  Report  R-100,  fejfca  corresponds  to  kif),  etc. 

The  integration  factor  for  airloads  is  defined  as 

'  2  X‘  IT  dri  ■  <45) 

SaR 


This  expression  is 
Report  R-100. 

Hie  integration 


the  same  as  that  given  by 


factor  for  dead-weight  loads 


k'sd 


=  2  So'  ~T^dr] 
Jo  V 


Eq.(F-l)  in 
is  de fined  as 
(46) 


This  factor  can  be  evaluated  from  a  knowledge  of  the  distri- 
butiort  of  the  dead  weight  of  wing  and  contents. 

Substituting  Eqs.(45)  and  (46)  into  Eq.(44),  the  volume 
of  shear  material  becomes 


(Vol),  “  W*ckitd)  • 


(47) 
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Factoring  out  the  terms  associated  with  airloads, 


(Vol), 


WgnfJs  bk isa  / 

*F°o  V 


"me 


(48) 


Equation  (48)  is  the  same  as  Eq.(36)  with  Jn  replaced  by  the 
term  in  parenthesis.  It  can  be  seen  that  this  term,  which 
represents  an  integrated  dead-weight  relief  factor,  reduces 
to  «/n(=l  —  Wwc/ Hg)  when  the  integration  factors  are  equal. 
This  would  occur  if,  as  previously  assumed,  the  dead  weight 
were  distributed  in  the  same  manner  as  the  airloads.  The  re¬ 
vised  inertia-load  integration  factor  for  shear  material 
will  be  denoted  by  Jnj ,  where 


j  _  i  _  k'sd 

6 


BENDING  MATERIAL 


The  volume  of  bending  material  which  is  proportional  to 
loading  is  given  by  Eq.(85)  of  Report  R-100: 

(Vol)fc  =  \  kibAbgb  (49) 


where  k^b  =  integration  factor  calculated  on  the  basis  of 
airloads 

A bg  =  cross-sectional  area  of  bending  material  re¬ 
quired  at  the  root  section,  and  where 


_  nf  (b/  2\  kb  Jn 

b*  "  ~  V’fcT/  Ua  ' 


(50) 


Substituting  Eq.(50)  into  Eq.(49),  the  volume  of  bending 
material  is  written 


(Vol)fc 


H'g'V  / b/2 \  bkb  kibJn 

6  \  /  keFa 


(51) 


Here,  as  in  Eq.(36),  Jn  gives  the  correct  value  for  dead¬ 
weight  load  relief  effects  only  at  the  root  section.  The 
span-wise  load  distribution  factor  kb  and  the  integration 
factor  kib  for  bending  material  are  based  only  on  airloads. 
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The  net  cross-sectional  area  of  bending  material  required 
at  any  station  is  given  by 

Ab  =  “  Md)  (52) 

ra  ke  h 

where  Ma  =  bending  moment  from  airloads 

Md  =  bending  moment  from  dead-weight  loads. 

The  volume  of  bending  material  is  found  by  integrating 
Eq. ( 52) : 

,  %  2b  ri  Ma  ,  2  6 

This  assumes  that  the  bending  moments  from  airloads  are 
greater  than  those  from  dead-weight  loads  at  all  span-wise 
sta  ti ons. 

The  value  of  Ma/h  at  the  root  section  is  given  by 

T*  *  7  *w(^K  (54> 


where  =  span-wise  load  distribution  factor  for  airloads. 

The  value  of  Md/h  at  the  root  section  is  given  by 

T1  =  7  (55) 


where  kf,d  =  span-wise  load  distribution  factor  for  dead¬ 
weight  loads. 

Multiplying  and  dividing  the  first  and  second  terms  of 
Eq .(53)  by  Eqs. (54)  and  (55),  respectively,  the  volume  of 
bending  material  is  written 

"wArV*.  ...  „ 


(voi)t  •  "vy  r  *, 

2  Fake  Jo  Ua/hR  1 


»«nf(^)bkbd 


,*e  J° 


i  Md/h 

todB/bR 


dr\  .  (56) 
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Since  the  allowable  axial  stress  is  constant,  the  cross- 
sectional  area  of  bending  material  required  at  any  station  is 
directly  proportional  to  the  quantity  M/h.  Therefore, 

j ~  =  (57) 

M°g/hR 

where  A(,o  =  area  of  bending  material  which  would  be  required 
to  carry  airloads  alone.  Also, 

4—  =  i T7T  <58» 

Kg  MdB/hR 


where  =  cross-sectional  area  of  bending  material  which 
would  be  required  to  carry  only  dead- weight  loads. 

Substituting  Eqs.(57)  and  (58)  into  Eq.(56)  and  factoring, 
the  volume  of  bending  material  becomes 


/  \  l/ty2\  brif  T  Ab  C i  Ab, 

(Vol)b  =  “(t  ) F  ,  ^gkba  I  ~  4  dr)  —  I  - 

/Fake  4'° 


dr,  I  .  (59) 


The  integration  factor  for  bending  material  required  by- 
airloads  is  defined  as 


^f'A 

*/o  Ak 


This  expression  is  the  same  as  the  original  equation  for  the 
bending-material  integration  factor  given  by  Eq.(F-6)  of 
Appendix  F,  Report  R- 100 . 

The  integration  factor  for  bending  material  which  would 
be  required  to  carry  only  dead-weight  loads  is  defined  as 


=  3  S'  7^ 

*/o  Afc 


Substituting  Eqs.(60)  and  (61)  into  Eq.(59),  the  volume 
of  bending  material  is  written 

(Vol)t  =  ("tKKa  ~  *.cKkiid) .  (62) 

/  ra  Re 
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Factoring  out  the  quantities  associated  only  with 
the  volume  of  bending  material  becomes 


ai rloads, 


(Vol)t 


tyg  nf  L  _  Ujcc  kbd  ^ijbd\ 

6  \fy?  )  Fa^e  \  Ifg  ^6a  iba ) 


(63) 


This  expression  for  volume  of  bending  material  is  the  same 
as  Eq.(5l),  except  that  the  dead-weight  relief  factor  J„  is 
replaced  by  the  integrated  relief  factor  in  parenthesis.  This 
revised  inertia-load  relief  factor  for  bending  material  is 
denoted  by  Jn,,  where 


/  _  i  _  ^ d 

1  T,  iT 


kikd 

kiba 


It  can  be  shown  algebraically  that  the  integration  factor 
for  shear  material  is  identical  with  the  span-wise  load 
distribution  factor  fej,.  For  purposes  of  explanation  the 
two  factors  are  left  in  their  original  form,  but  this 
identity  should  be  kept  in  mind  when  calculating  actual 
values. 


Some  actual  values  for  the  integration  factors  and  relief 
factors  are  given  below: 


Airplane 

_ 

Item 

A 

B 

!■  *i.a 

.810 

.806 

2>  *‘sd 

.599 

.565 

3-  fc»b« 

1.  128 

1.030 

*‘bd 

.706 

.813 

5.  J„ 

.  546 

.647 

6-  Jnm 

.665 

.752 

7  •  Jnb 

.790 

.805 

8.  Percent  increase  in  weight  of 

shear  material  by  use  of 
corrected  relief  fsctor 

22 

16 

9.  Percent  increase  in  weight  of 
bending  meterial  by  use  of 
corrected  relief  factor 

45 
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Designation  of  airplanes: 


A  =  largo  four-engine  cargo  airplane  calculated  for  the 
critical  condition  of  no  wing  fuel 

B  =  fighter  aircraft  aith  noenginaain  the  eiaga  calcu¬ 
lated  for  tha  critical  coaditioa  of  no  wing  fuel. 
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The  errors  introduced  by  using  Jn  rather  than  th 
grated  relief  factor  can  be  seen  by  comparing  Items 
with  Item  5.  The  percentage  error  is  seen  to  be  lar 
the  four-engine  cargo  airplane. 


e  inte- 
6  and  7 
ger  for 


When  more  detailed  data  are 
actual  airplanes,  these  factors 
er  of  airplanes,  and 
d  for  different  clas 


for  a  numb 
be  compute 


acquired  on 

(Jn,  and  Jnb) 
perhaps  ave 
ses  of  airpl 


wing  loads  for 
can  be  evaluated 
rage  values  can 
anes . 
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APPENDIX  B 

DERIVATION  OF  THE  INTEGRATED  SHEAR 
REDUCTION  FACTOR 


The  expression  for  the  shear  reduction  factor  Js  given 
byEq.(A-8)of  Report  R-100  is  valid  only  at  the  root  section 
of  the  wing.  An  investigation  of  actual  airplane  wings  has 
shown  that  the  variation  in  the  value  of  Js  along  the  span 
makes  necessary  the  use  of  an  average  or  "integrated"  value. 
This  integrated  shear  reduction  factor,  when  substituted 
in  the  general  wing  weight  equation,  yields  an  amount  of 
shear  material  which  is  corrected  for  the  over-all  effect 
of  shear  relief  along  the  span. 

The  integrated  shear  reduction  factor,  denoted  by  .  is 
defined  as  follows: 


_  ( Vo  1  )Sc 

'  (Voi)Sn 


where  (Vol  )s<_  =  volume  of  shear  material  necessary  to  carry 

the  shear  load  after  it  is  corrected  for  re¬ 
lieving  shear  loads  from  the  flanges 

(Vol)jQ  =  volume  of  shear  material  necessary  to  carry 
the  total  uncorrected  shear  load. 

This  expression  may  also  be  written 


j;  =  i  - 


where  (Vol),^  =  volume  of  material  which  is  subtracted  from 

so  when  relieving  shear  is  considered. 

The  volume  of  shear  material  before  correcting  for  shear 
relief  is 


(Vol  )*0  =  ~ 

0  9  F. 

*  *o 


where  VK  =  maximum  shear  load  at  the  root  of  one  wing. 
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Substituting  Eq.(73)  into  £q.(67)  and  expressing  the 
span-wise  distance  y  in  terms  of  7),  where  r\  =  y/(b/2), 


(Vol)s,  = 


2(1  ~  «)  ri  _ 

•'o  1  - 


T)  +  Tjm 


Multiplying  and  dividing  by  the  bending  moment  A/fl  at  the 
root  section, 


/„  ,n  2d-*)%  r1  A 

(Vol)S/  =  — — -  I  — 

1  keFSn  j  0  1- 


r)  +  T}m 


Equations  (F-5)  and  (F-6)  of  Appendix  F,  Report  R-100, 
show  that 


P- 

Jo  1  - 


M/Mr  fci b_ 

-  7]  +  r\m  3 


where  kih  =  integration  factor  for  bending  material. 
Substituting  Eq.(76)  into  Eq.(75), 


(Vol)4 


__  2  (1  ~m)MRklb 


Substituting  Eqs. (77 )  and  (66 )  into  Eq. ( 65) ,  the  integrated 
shear  reduction  factor  becomes 


-  l  —  4/%\  £j6(l 

S  3 \vj  kekls 


As  shown  in  Appendix  F  of  Report  R-100,  the  following 
expression  may  be  written 


=  *6  7  =  *<.7 


Substituting  Eq.(79)  into  Eq.(78)  and  simplifying,  the 
integrated  shear  reduction  factor  is  written 

Js  =  1  ~  k,t(1~m)  .  (80) 

3ke 


The  expression  for  J,'  given  by  Eq.(80)  is  derived  here 
to  replace  Eq.(A-8)  of  Report  R-100.  Equation  (80)  has  been 
derived  on  the  basis  of  net  shear  and  bending  loads.  When 
the  airloads  and  inertia  loads  are  treated  separately,  it 
is  necessary  to  go  through  the  derivation  above  and  to 
express  the  net  moment  M  as  Ma  —  ,  the  shear  V  as  Va  —  Frf, 
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etc.  (The  subscript  a  denotes  airload  values  and  d  denotes 
dead-weight  values. )  Such  a  derivation  gives  the  expression 
for  Jj  as 


Js 


1  *t6„(l  ~  m> 


3fee 


Jnb 


(81) 
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APPENDIX  C 


DERIVATION  OF  REVISED  RIB  WEIGHT  EXPRESSIONS 


The  original  expressions  for  the  volume  of  rib  material 
are  given  in  Section  II  as  Eqs.(le)  through  (lg).  The 
derivations  of  these  expressions  are  given  in  Report  R-100. 
The  derivation  for  the  new  rib  weight  terms  is  given  below. 

RIB  SHEAR  MATERIAL 

The  expression  for  the  area  (average  volume  per  foot 
of  span)  of  rib  shear  material  is  given  by  Eq.(47)  of  Report 
R-100  as 

£*2  +  iai3lK  (82) 

L  2  (H4S)F,„ 

Since  the  term  Cs  is  the  intercept  (at  V/h2  =  0)  value  of 
t/h  (see  Fig. 5,  Report  R-100),  it  represents  a  minimum  al¬ 
lowable  ratio  of  average  thickness  to  depth.  A  study  of 
actual  airplanes  showed  that  it  is  desirable  to  place  the 
minimum- value  limitation  on  thickness  rather  than  on  proportion. 

Since  Cs  =  ( t/h)0 ,  the  minimum  average  thickness  lrQ  for 
the  rib  shear  material  is  given  by 

tr  =  Csh  .  (83) 

ro  a 

Substituting  Eq.(83)  into  the  first  term  of  Eq.(82),  the 
expression  for  the  average  area  (volume  per  unit  length  of 
span)  of  rib  shear  material  becomes 


a'  - 1  hB  + i  nfJn 

rs  r°  L  2  ( 144.S  )fs 


The  volume  is  found  by  integrating  Eq. (84)  over  the  span, 
as  in  Report  R-100: 

(Vo l)r,  =  JrnKhKB(H_4)S  +  1  Jnnf*gKBC  .  {g5) 

'  4  1 /  rtf 
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RIB  FLANGE  MATERIAL 


As  previously  stated,  the  new  expression  for  the  volume 
of  rib  flange  material  is  derived  on  the  basis  of  axial 
loads  associated  with  chordwise  bending  moments. 

From  Fig. 10  and  Eq.(45)  of  Report  R-100,  it  can  be  seen 
that  the  expression  for  the  maximum  chordwise  bending 
moment  is  given  by 

1  V 

^max  —  ~  nf  J n  '  L  CB  •  (86) 

o  1 44o 


The  axial  load  P  in  the  rib  flanges  is  given  by 

_  M  _  nfJn  WgLCB 
~h  8  h(  144S) 


(87) 


The  cross-sectional  area,  per  bay,  of  rib  flange  material  is 


2P  _  nfJn  WgLCB 


Fy  4h(144S)Fy 


(88) 


where  Fy  =  compressive  yield  stress,  or  equivalent  allowable 
stress . 


The  volume  of  rib  flange  material  per  bay  is 
( Vo  1 )  r  /  =  AB  =  nfJnWgLCB 2 


bay 


4/i(  1445  )FV 


(89) 


Expressing  L,  B,  and  h  as  fractions  of  chord,  the  average 
area  (average  volume  per  unit  length  of  span)  is  written 


A'rf  = 


_  nfJ„  WgR^C4 


4rKh(lUS)Fy 


(90) 


Integrating  Eq. (90)  over  the  span  gives  the  total  volume 
of  rib  flange  material  as 


(Vol)r/ 


nfJn  WgKjC 
4  KhFy 


(91) 


34 


APPENDIX  D 


DERIVATION  OF  EXPRESSIONS  FOR  COMBINING 
INTEGRATION  FACTORS 


Since  integration  factors  are  calculated  separately  for 
various  concentrated  and  distributed  loads,  they  must  be 
combined  into  single  factors  in  order  to  be  used  in  the 
wing  weight  equations.  This  may  be  done  by  means  of  the 
simple  expressions  derived  below. 


SHEAR  INTEGRATION  FACTORS 

The  shear  integration  factor  represents,  in  nondimensional 
form,  the  area  under  the  shear  curve.  Figure  6  shows  the 
shear  curves  for  two  sources  of  load,  as  well  as  the  total 
shear. 


Fig.  6 — Shear  diagram  for  two  superimposed 
loading  conditions 

The  integration  factor  corresponding  to  shear  curve 
V,  is  given  by 

l  =  - il_  (92) 


1,1  1  /6\ 
iMi) 


where  Aj  =  area  under  curve  of  Vj . 

Similarly,  the  integration  factor  for  is  written 


*  e  1 


where  A]  =  area  under  curve  of  V% . 
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The  over-all  integration  factor  ^isT,  which  is  to  be 
written  in  terms  of  and  ^isa*  is  given  by 


kiT  =  - 


A\  +  4  2 


o  (K*a  +If*.) 


(l) 


(94) 


Solving  Eqs.(92)  and  (93)  for  <4j  and  A  2,  respectively,  and 
substituting  into  Eq.(94),  the  total  integration  factor 
becomes 


i*r 


(95) 


Simplifying,  the  equation  may  be  written  in  general  form  as 


+ 


+  F. 


(96) 


This  derivation  assumes  that  the  total  shear  does  not  change 
sign  at  any  point  along  the  span. 


BENDING  INTEGRATION  FACTORS 

The  bending  integration  factor  represents,  in  non- 
dimensional  form,  the  area  under  the  curve  of  M/h,  which 
is  the  flange  load.  Figure  7  shows  curves  of  M/h  for  two 
sources  of  bending  loads  and  a  curve  for  the  total  value 
of  M/h. 
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Fig.  7 — Diagram  of  axial  flange  loads  for  two  sources 
of  bending  moments 


This  derivation  assumes  that  the  total  bending  moment  does 
not  change  sign  at  any  point  along  the  span. 
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